The a — > 7 transition in Ce: a theoretical view from optical spectroscopy. 
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Using a novel approach to calculate optical properties of strongly correlated systems, we address 
the old question of the physical origin of the a — > 7 transition in Ce. We find that the Kondo collapse 
model describes the optical data better than the Mott transition picture. Our results compare well 
with existing thin-film experiments. We predict the full temperature dependence of the optical 
spectra and find the development of a pseudogap around 0.6 eV in the vicinity of the a — > 7 phase 
transition. 

PACS numbers: 71.27.+a,71.30.+h 



At a temperature of 600 K and pressure less than 
20 Kbars, elemental Cerium undergoes a transition be- 
tween two isostructural phases, a high pressure phase or 
a phase and a low pressure 7 phase. In a—Ce the / elec- 
tron is delocalized (for example the spin susceptibility 
is temperature independent ) while in 7-Ce the / elec- 
tron is localized (for example the electron has a Curie- 
like susceptibility). Several basic questions about this 
transition are still being debated. What is the driving 
mechanism of this transition? What is the role of the 
spd electrons which are near the Fermi level in this ma- 
terial? To address these questions two main hypothesis 
have been advanced. Johansson proposed a Mott tran- 
sition scenario Q, where the transition is connected to 
delocalization of the / electron. In the a phase the / 
electron is itinerant, band like while in the 7 phase it 
is localized and hence does not participate in the bond- 
ing, explaining the volume collapse. In this picture the 
spd electrons are mere spectators well described by the 
density functional theory calculations (DFT) in the local 
density approximation (LDA) or its extensions. The / 
electrons in a-Ce can also be treated in the band picture 
while in 7-Cerium should be considered as part of the 
core, and do not participate in the bonding. 

A different view on this problem was proposed by Allen 
and Martin who introduced the Kondo volume collapse 
model for the a — 7 transition. They suggested that the 
transition was connected with modifications in the effec- 
tive hybridization of the spd bands with the /-electron. 
In this picture the main change when going from a to 
7 is the degree of hybridization and hence the Kondo 
scale. In a series of publications Q they implemented 
this idea mathematically by estimating free energy dif- 
ferences between these phases by using the solution of 
a Anderson-Kondo impurity model supplemented with 
elastic energy terms. This approach, views the lattice 
of Cerium atoms in terms of a single impurity Ander- 
son model, describing the Ce / electron in fixed bath 
of band electrons representing the spd electrons. The 



modern Dynamical Mean Field Theory (DMFT) mod- 
ifies this simple single impurity picture by imposing a 
self-consistency condition on the bath that the / elec- 
trons experience. The full power of the DMFT method 
in combination with realistic band structure calculations 
(LDA+DMFT) was recently brought to bear on this 
problem LDA+DMFT allowed the computation of 

the photoemission spectra of Cerium in both phases, and 
the thermodynamics of the transition starting from first 
principles. The theoretical results were in good agree- 
ment with existing experiments |fj. The photoemission 
spectra close to the Fermi level, is dominated by the / 
electron density of states. The most recent DMFT stud- 
ies of Cerium demonstrated that both the Mott transition 
picture and the Kondo collapse picture lead to similar 
photoemission and inverse photoemission spectra. The 
spectra of the a phase consists of Hubbard bands and 
a quasiparticle peak while the insulating 7 phase has no 
quasiparticle or Kondo peak in the spectra and consists 
of Hubbard bands only. Hence, from the study of the 
/ electron spectra which dominates the photoemission 
spectroscopy, it is not possible to decide whether the res- 
onance arises from the proximity to the Mott transition 
or from the Kondo effect between the / and the the spd 
electrons. 

In this letter we revisit this problem theoretically using 
optical spectroscopy. Our qualitative idea is that the spd 
electrons have very large velocities, and therefore they 
will dominate the optical spectrum of this material. In 
the Mott transition picture, the spd electrons are pure 
spectators, and hence no appreciable changes in the op- 
tical spectrum should be observed. On the other hand, 
if the hybridization between the spd electrons and the / 
electrons increases upon entering the a phase, we expect 
a pseudogap to develop as the temperature is lowered 
because spd carriers are strongly modified as they bind 
to the / electrons. Optical experiments on thin films 
were carried out recently Q . They offer a benchmark to 
compare the results of the theoretical DMFT calculation. 
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We will show that DMFT reproduces all the features ob- 
served in the experiment. By studying in detail the tem- 
perature dependence of the spectra, we can interpret the 
different experimental features. We find that Cerium is 
much closer to the volume collapse picture than to the 
Mott transition picture. 

Formalism - Within the LDA+DMFT method 0, the 
LDA Hamiltonian is superposed by a Hubbard-like lo- 
cal Coulomb interaction, which is the most important 
source of strong correlations in correlated materials and 
is not adequately treated within LDA alone. The result- 
ing many-body problem is then treated in DMFT spirit, 
i.e., neglecting the nonlocal part of self-energy. It is well 
understood by now, that this theory is not only exact in 
the limit of infinite dimensions but is also a very valuable 
approximation for many three dimensional systems since 
it is capable of treating delocalized electrons from LDA 
bands as well as localized electrons on equal footing. 

To calculate the central object of DMFT, the local 
Green's function Gi oc , we solved the Dyson equation 

{Ht DA + E(«) - E dc - enMlg* = 0, (1) 

where H^ DA is the LDA Hamiltonian expressed in a lo- 
calized linear muffin-tin orbital (LMTO) base, Ok is the 
overlap matrix appearing due to nonorthogonality of the 
base and is the self-energy matrix to be determined 
by the DMFT. A double counting term Edc appears here 
since the Coulomb interaction is also treated by LDA in 
a static way therefore the LDA local correlation energy 
has to be subtracted. Given the eigenvalues €kj U , the 
left eigenvectors ip k ^ and the right eigenvectors ?/> k j/ of 
Eq. the local Green's function can be expressed by 
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The local self-energy E(w), that appears in Eq. (|TJ, can 
be calculated from the corresponding impurity problem, 
defined by the DMFT self-consistency condition 



Gioc — (cj — E imp — £(u>) — A imp (co)) 



(3) 



Here A imp is the impurity hybridization matrix and E. 
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purity problem, i.e., the functional E[Aj mp (o;), Ei mp , U], 
closes the set of equations © an d ©■ 

Various many-body techniques can be used to solve 
the impurity problem, among others, the Quantum 
Monte Carlo method (QMC), Non-crossing approxi- 
mation (NCA) or iterative perturbation theory (IPT). 
Here, we used One-crossing approximation (OCA) [3, 
licj . which is, like NCA, a self-consistent diagrammatic 
method that perturbes on atomic limit. Unlike NCA, 
it takes into account certain type of vertex correction, 
namely all diagrams that do not contain a line with more 
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FIG. 1: Luttinger-Ward Functional for the impurity solver 
called one-crossing approximation (OCA). Full lines corre- 
spond to the conduction electrons from the bath, while dashed 
lines stand for slave particles. Each slave particles, however, 
corresponds to exactly one atomic states. 



than one crossing, hence the name, One crossing approxi- 
mation. This approximation is far superior to NCA, as is 
well known, considerably improves the Fermi-liquid scale 
Tk and reduces some pathologies of NCA, but is more 
time consuming. The OCA is also the lowest order self- 
consistent approximation exact up to V 2 oc Aj mp where 



are the impurity levels. The solution of the Anderson im- 



is the impurity hybridization with the electronic 
bath. A self-consistent diagrammatic approximation is 
most conveniently expressed by its Luttinger-Ward func- 
tional $, which uniquely determines the dressing of the 
slave particles, i.e., S Q/ g = d&/dGp a . First diagram of 
<I> in Fig. 2] corresponds to the well known NCA, while 
the second, i.e., the crossing diagram, defines corrections 
that are taken into account in OCA. 

The optical conductivity is given by a current-current 
correlation function which can be expressed in the fol- 
lowing way 
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where we have denoted = e±u;/2, and used the short- ments M^jf appear as standard dipole allowed transition 
cut notations ej. g = £kje, e kj£ = € t.js- ^ ne ma "trix ele- probabilities which are now defined with the right and left 
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solutions ip R and ip L of the Dyson equation: 

CK1 ,0!2 

£ (Vi' £ < Q 3) S '<a3a 4 (^--,« 4 ) S >(5) 



where we have denoted V'ke'" 



and 



assumed that (V'ke) + = 



while (v£t = (iter- 

The expressions |@J,l(5J represent generalizations of the 
formulae for optical conductivity for a strongly correlated 
system, and involve the extra internal frequency integral 
appearing in Eq. 

As in previous work 0, 0| we treat only the /-electrons 
as strongly correlated thus requiring full energy resolu- 
tion, while all other electrons such as Ce s,p,d are as- 
sumed to be well described by the LDA. The spin-orbit 
coupling is fully taken into account therefore the size of 
the matrixes in Eq. (JJJ is 32 x 32, while the self-energy 
matrix appears in a sub-block of size 14 x 14. The value 
of the Coulomb interaction U was calculated by the lo- 
cal density constrained occupation calculations 01 an( i 
found to be U ~ 6 eV. After the self-consistency of the 
dynamical mean field equations is reached, the quasipar- 
ticle spectra described by i/^ and ek JW are found, and 
the expression 10} for cf^^uj) can be evaluated. Here, 
we paid a special attention to the energy denominator 
l/(u> + ejC _ — e kj'e+ ) appearing in Due to its strong 
k-dependence the tetrahedron method of Lambin and 
Vigneron is used. The integral over internal energy 
e is calculated in an analogous way. The frequency axis is 
divided into discrete set of points and assumed that the 

eigenvalues ekj e and the matrix elements ', >iV [s^ , e + ) 
can be linearly interpolated between each pair of points. 
We first perform the frequency integration to convert the 
double pole expression (@J into a single pole expression 
for which the tetrahedron method is best suited. 

Results - What distinguishes the a phase from the 7 
phase when looking at the optical response of Cerium? 
How do the calculated results compare with experiments? 
In the top panel of Fig. |2 we present the calculation of 
LDA+DMFT method while in the bottom panel the re- 
sults of experiment on thin films Q are reproduced. No- 
tice that in both the experiment and theory the overall 
magnitude of conductivity in 7 phase is smaller than in 
a due to the larger volume and hence reduced velocity 
of 7 phase. In the a phase, the Drude peak width is 
of the order of 0.1 eV while in 7 phase it is at least an 
order of magnitude narrower. This provides us with a 
crude estimate of Tk, which is of the order of 1000 K 
and 100 K in the a and 7 phase, respectively. The inter- 
mediate frequency range in a is characterized by a dip 
or pseudogap ranging from 0.2 eV to 0.8 eV and a peak 
around 1 eV. This gap is filled-in in 7 phase resulting 
in a broad hump up to 0.5 eV. This nontrivial behavior 
can be understood by examining the orbitally resolved 
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FIG. 2: The top panel shows the optical conductivity cal- 
culated by LDA+DMFT method for both a and 7 phase of 
Cerium. The volume of a phase is 28.06A 3 and the temper- 
ature used in calculation is 580K. The volume of 7 phase is 
34.37A and T = 1160-fC. The bottom panel shows experi- 
mental results measured by the group of D. van der Marel 0. 
The measurements for a phase were done at 5 K and for 7 
phase at 300 K. 
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FIG. 3: Temperature dependence of partial density of states 
for both phases a and 7 of Ce. 



partial density of states displayed in Fig. 

The largest contribution to the total density of states 
around the Fermi level comes from the / bands which 
have much smaller velocities and, at the same time, al- 
most all spectral weight above the Fermi level. Hence, 
the transitions from below to above the Fermi level have 
a small amplitude therefore / bands make a very small 
contribution to the optical response. Indeed, we found 
that the contribution involving only / electrons is orders 
of magnitude smaller than contribution of the d bands. 
Since s and especially p bands have very little spectral 
weight close to the Fermi level, the most important con- 
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FIG. 4: The partial contribution to the optical conductivity 
of a phase of Ce at T = 11QQK. 



tribution to the optical response comes from the d bands. 

To shed some light on the states probed optically we 
introduce an orbital resolution of the optical conductivity 
(or "fat optics"). 
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where // stands for the pure / contribution, dd for pure 
spd contribution and fd for the mixed part. The partial 
conductivities are calculated by taking the appropriate 
dipole matrix elements M^jji in Eq. ifljl. for example, to 
calculate // part all four indexes a>i, a.%, 0:3, a4 in Eq. (J5J 
run over / orbitals only while in dd part all indexes run 
over spd orbitals. The results displayed in Fig.0J demon- 
strate that indeed the optical conductivity has very little 
/ character. 

The small contribution of the / electrons to the opti- 
cal conductivity, however, does not imply that the large 
Kondo peak, which has / character, is irrelevant for the 
optical conductivity. The / electrons, mixing with the 
spd bands, form a large Kondo peak that increases as 
the temperature is reduced and eventually induces a hy- 
bridization pseudogap in spd spectra. As one can see in 
Fig. |3 the d bands in a phase have a very pronounced 
pseudogap which is growing by lowering temperature ex- 
actly as the Kondo peak builds up. The spectral weight 
is transferred from the Fermi level into the side-peaks 
which are 1 eV apart causing 1 eV peak in optical con- 
ductivity. In the 7 phase, the Kondo peak disappears 
because the effective hybridization of the / with spd elec- 
trons is smaller and Kondo scale Tk is reduced for at least 
an order of magnitude. In this regime the hybridization 
pseudogap in the d bands is not formed, and instead we 
find a broad hump in the optical response. 

Although all basic features of measured optical re- 
sponse are well explained by LDA+DMFT calculation, 
there are some discrepancies. The overall magnitude of 
measured conductivity in both phases is approximately 
half of the calculated value. This discrepancy might be 
a consequence of negligence of interactions on the spd 
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FIG. 5: Temperature dependence of optical conductivity for 
both phases a and 7 of Ce. 



electrons, an effect which would reduce the current ma- 
trix elements. Early measurements by Rhee. et.al. jg], 
hovewer, suggest that the optical conductivity is indeed 
for factor of two larger than the one measured by the 
group of D. van der Marel [3. Also the shoulder in a 
phase, appearing around 0.33 eV in Fig. 03, is absent in 
LDA+DMFT results. 

The temperature dependence of the optical response 
for both the a and 7 phase are shown in Fig. |3J With 
increasing temperature, the 1 eV peak moves to smaller 
frequencies and finally at few 1000 K disappears. At the 
same time, the pseudogap below the peak gradually dis- 
appears as the temperature is raised and evolves into a 
broad hump, just like the one in 7 phase. This temper- 
ature dependence is easily understood by looking at the 
partial density of states in Fig. [3] Since the Kondo peak 
is gradually reduced with increasing temperature the hy- 
bridization pseudogap in d bands disappears causing fea- 
tureless optical response. In 7 phase, the Kondo peak is 
very small and causes some temperature dependence at 
very low temperatures only at very low frequency below 
0.1 eV. 

Conclusion - The methodology we introduced, allow us 
to interpret the optical spectra of a and 7 Cerium, in fa- 
vor of the Kondo volume collapse model. The formalism 
should be useful in many other problems of correlated 
electrons where temperature dependent transfer of spec- 
tral weight between widely different energy ranges takes 
place, a phenomena which can be beyond the scope of 
simple model Hamiltonians. 



[1] B. Johansson, Philos. Mag. 30, 469 (1974). 

[2] J.W. Allen and R.M. Martin, Phys. Rev. Lett. 49, 1106 



5 



(1982); J.W. Allen and L.Z. Liu, Phys. Rev. B 46, 5047 
(1992). 

[3] For recent reviews see K. Held et al, Psi-k Newslet- 
ter #56 (April 2003), p. 65; A. I. Lichtenstein, M. 
I. Katsnelson, and G. Kotliar,in Electron Correlations 
and Materials Properties 2, ed. A. Gonis (Kluwer, 
NY) cond-mat/02 11076 . 

[4] M.B. Zolfi, 1 A. Nekrasov , Th. Pruschke, V.I. Anisimov, 
and J. Keller, Phys. Rev. Lett 87, 276403 (2001). 

[5] K. Held, A.K. McMahan, and R.T. Scalettar, Phys. Rev. 
Lett. 87, 276404 (2001); A.K. McMahan, K. Held, and 
R.T. Scalettar, Phys. Rev. B 67, 075108 (2003). 

[6] D.M. Wieliczka, C.G. Olson, and D.W. Lynch, Phys. 
Rev. Lett. 52, 2180 (1984); E. Weschke, C. Laubschat, T. 
Simmons, M. Domke, O. Strebel, and G. Kaindl, Phys. 



Rev. B 44, 8304 (1991); F. Patthey, B. Delley, W.D. 

Schneider, and Y. Baer, Phys. Rev. Lett. 55, 1518 (1985). 
[7] J.W. van der Eb, A. B. Kuz'menko, and D. van der Marel, 

Phys. Rev. Lett. 86, 3407 (2001). 
[8] J.Y. Rhee, X. Wang, B.N. Harmon, and D.W. Lynch, 

Phys. Rev. B 51, 17390 (1995). 
[9] Th. Pruschke and N. Grewe, Z. Phys. B: Condens. Matter 

74, 439, 1989. 

[10] K. Haule, S. Kirchner, J. Kroha, and P. Wolfle, Phys. 

Rev. B 64, 155111, (2001). 
[11] Ph. Lambin and J. P. Vigneron, Phys. Rev. B 29, 3430, 

(1984). 

[12] A.K. McMahan, C. Huscroft, R.T. Scalettar, and E.L. 
Pollock, J. Comput.-Aided Mater. Des. 5, 131 (1998). 



